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This work establishes a firm relationship between classical nonlinear resonances and the phe- 
nomenon of dynamical tunneling. It is shown that the classical phase space with its hierarchy of 
resonance islands completely characterizes dynamical tunneling and explicit forms of the dynamical 
barriers can be obtained only by identifying the key resonances. Relationship between the phase 
space viewpoint and the quantum mechanical superexchange approach is discussed in near-integrable 
and mixed regular-chaotic situations. For near-integrable systems with sufficient anharmonicity the 
effect of multiple resonances i.e., resonance-assisted tunneling can be incorporated approximately. 
It is also argued that the, presumed, relation of avoided crossings to nonlinear resonances does not 
have to be invoked in order to understand dynamical tunneling. For molecules with low density of 
states the resonance-assisted mechanism is expected to be dominant. 
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I. INTRODUCTION 

Intramolecular vibrational energy redistribution (IVR) 
in a molecule, from an initially prepared nostationary 
state, is at the heart of chemical reaction dynamics. The 
timescales and mechanisms involved in this process of 
energy flow have been investigated in great detail by ex- 
periments and theory^. Although considerable progress 
has been made over the years there are many aspects 
of this phenomenon whose understanding still eludes us. 
Nevertheless it is expected that there is a part of the 
IVR mechanism which can be usefully understood based 
on classical dynamics alone and there is a part of IVR 
which is intrinsically quantum mechanical. Both classical 
and quantum mechanisms coexist in a given molecule and 
give rise to complicated spectral patterns and splittings. 

It is now well established that a specific overtone ex- 
citation in a molecule will undergo IVR if the initial 
nonstationary state, zeroth-order bright state (ZOBS), 
is coupled to other 'dark' zeroth-order states via strong 
anharmonic resonances. In such cases classical dynam- 
ics, at various levels of detail and sophistication, can pro- 
vide useful insights into the timescales and mechanisms of 
IVR^. On the other hand if the ZOBS is not close to any 
resonance then classically one expects no IVR and hence 
no fractionation of spectral lines. However, it is still pos- 
sible for IVR to occur via quantum routes. In other words 
the initial state will mix with other states giving rise to 
spectral splittings and, possibly, complicated eigenstates. 
Since the mixing is classically 'forbidden' it would be ap- 
propriate to associate some kind of tunneling with such 
quantum routes to energy flow. Indeed such a suggestion 
was made over two decades agc'^ and the term 'dynam- 
ical tunneling' was coined. The notion of tunneling is 
meaningless without one form of a barrier or other and 
the term 'dynamical' is prefixed to distinguish from the 
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usual coordinate space tunneling through static potential 
barriers. The barriers in dynamical tunneling, in general, 
are more subtle to identify and exist in the phase space. 

Dynamical tunneling can have important consequences 
for the interpretation of molecular spectra since the fin- 
gerprints of IVR are spectrally encoded in the form of 
intensities and splittings. Traditionally dynamical tun- 
neling in the molecular context has been associated with 
the occurence of local mod^ stretches in symmetric 
molecules. However it is important to emphasize that 
the concept of dynamical tunneling is more general - 
any flow of quantum probability between regions which 
are classically disconnected is dynamical tunneling. The 
identification of classically disconnected regions requires 
the knowledge of the phase space topology and hence 
it is not very surprising that dynamical tunneling is in- 
evitably linked to the underlying classical dynamics. In 
the context of local mode doublets in symmetric systems 
important work by Jaffe and Brumer^ and Kellman'' pro- 
vided classical phase space perspectives on the normal to 
local transition. One of the reasons for the interest in the 
near-degenerate local doublets has to do with their long 
lifetimes. Thus provided one can experimentally prepare 
a high overtone, of say the OH-stretch in H2O, then the 
extremely small splitting provides a large window of time 
to perform mode selective chemistry. The main prob- 
lem with such an idea is that at such high vibrational 
excitations other Fermi resonances and rovibrational in- 
teractions could destroy the degeneracies^. We hope to 
shed some light on these issues in this work wherein the 
system of interest does exhibit close degeneracies despite 
multiple resonances and a classical phase space which is 
mixed regular-chaotic. Interestingly, similar issues have 
been addressed in the context of the existence of discrete 
breathers in a network of nonlinear oscillators^. 

One of the earliest examples in symmetric systems are 
the local mode doublets observed in the water molecule 
which were explained by Lawton and Child as due to dy- 
namical tunneling^. The experiments^*' of Kerstel et.al. 
on unusally slow (hundreds of picoseconds) intramolec- 
ular vibrational relaxation of CH-stretch excitation in 
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(CH3)3CCCH molecule has also been ascribed to dy- 
namical tunneling by Stuchebrukhov and Marcusii. In 
this case the extremely slow IVR out of the CH-stretch 
is an instance of dynamical tunnneling in nonsymmet- 
ric systems. In a later experiment Mcllroy and Ncs- 
bitt noted multiquantum state mixing due to very small 
couplings (0.1 cm"-'^ or less) in the CH stretch excita- 
tion in propynei^. Recent experiments have shown sim- 
ilar effects in the CH-stretch excitations of triazinei^ 
and 1-butyne^- and there are reasons to expect dy- 
namical tunneling to manifest in general rovibrational 
spectrai^Oeiios^^ 

Much of the theoretical understanding of dynamical 
tunneling has come about from the analysis of bent tri- 
atomic ABA molecules. The focus primarily has been on 
vibrations and this work is no expception in this regard. 
However we note the important work by LehmannAS 
wherein the nontrivial effects arising due to coupling be- 
tween local modes and rotations are studied in detail. 
Davis and Heller'^ emphasized a phase space picture and 
implicated classical resonances in the phase space as the 
agents of dynamical tunneling. In this approach classical 
trajectories trapped in one region of the phase space were 
imagined to be separated by dynamical barriers, due to 
the resonances, from symmetrically equivalent regions of 
the phase space. The degeneracy is then broken by dy- 
namical tunneling. A clear demonstration of the role of 
isolated resonances was provied by Ozorio de Almeida as 
well2£. A slightly different picture was provided by Sib- 
ert, Reinhardt and Hynes^^ in their work on energy flow 
and local mode splittings in the water molecule. The set- 
ting was again in terms of classical nonlinear resonances 
but the dynamical barrier was identified in angle space 
as opposed to in action space. Later Hutchinson, Sib- 
ert and Hynes-- provided an explanation for the quan- 
tum energy flow in terms of high order perturbation the- 
ory. Stuchebrukhov and Marcus^i reanalyzed the ABA 
system in terms of chains of off-resonance virtual states 
(" vibrational superexchange" ) connecting any two degen- 
erate local mode states. An important result was the 
equivalence between the vibrational superexchange ap- 
proach and the usual WKB expression for splittings in 
a doublewell system. This equivalence between pertur- 
bative and the WKB approaches was shown for a classi- 
cally integrable case. Demonstration of such an equiva- 
lence in more general situations is an open problem. De- 
spite seemingly different approaches a crucial ingredient 
to any description of dynamical tunneling involves the 
various anharmonic resonances. Depending on the level 
of excitation the underlying phase space can exhibit sev- 
eral (near) equivalent regions separated by nonlinear res- 
onances (near-integrable phase space) or chaos (mixed 
phase space). For large molecules with sufficient den- 
sity of states, in the near-integrable regime. Heller has 
conjectured-^ that a nominal 10~^-10~^ cm""'^ broaden- 
ing of spectroscopically prepared zeroth-order states is 
due to dynamical tunneling between remote regions of 
phase space facilitated by distant resonances. 



The aforementioned conjecture is based on the notion 
that phase space is the correct setting for an understand- 
ing of dynamical tunneling. It is thus natural to expect 
the splittings to be sensitive to the structure of the phase 
space. It also follows that given the phase space structure 
one ought to be able to compute the associated splittings. 
General forms for splittings cannot be written down eas- 
ily since explicit forms of the dynamical barriers in the 
phase space, seperating quantum states localized in dis- 
tant regions of the phase space, can only be provided 
upon identification of the key nonlinear resonances. This 
is also tanatamount to identifying the mechanism of dy- 
namical tunneling and hence IVR. Two factors make this 
a difficult problem. Firstly, our understanding of global 
structure of multidimensional classical phase space is still 
in its infancy. Thus without a global view of the phase 
space at energies corresponding to the doublets it is dif- 
ficult to identify the main nonlinear resonances and per- 
haps the presence of appreciable stochastic layers. Sec- 
ondly, and related to the first factor, important work over 
the last decade has established that dynamical tunnel- 
ing is sensitive not only to the nonlinear resonance o^^i^^ 
but also to the classical stochasticity^^. In the near- 
integrable regimes dynamical tunneling can be enhanced 
by many orders of magnitude due to the various non- 
linear resonances. This is called as resonance-assisted 
tunneling and it has been suggested^^ that the nonlinear 
resonances play an important role in mixed phase space 
scenario as well. In cases where an appreciable chaotic 
region separates the two symmetry related regular zones 
one has the phenomenon of chaos-assisted tunneling. The 
hallmark of chaos-assisted tunneling is the erratic fluctu- 
ations of the splittings^- due to 'chaotic' states having 
nonzero overlaps with the regular doublets. In particular 
the splittings show algebraic dependence on h as opposed 
to the integrable exp(— l/7i) scaling. Experiments^ on a 
wide variety of systems have highlighted the role of the 
phase space structures in dynamical tunneling. 

Is the evidence for the sensitivity of dynamical tunnel- 
ing to the phase space structures already present in the 
molecular spectra? The answer to this question, apart 
from a fundamental viewpoint, is also potentially rel- 
evant to mode-specific chemistry and control. In this 
context, it is significant to note that most of the earlier 
works in the molecular context have focused on dynami- 
cal tunneling where only one resonance was involved with 
no chaos. A notable exception is the work of Davis and 
Heller— wherein a hint to the role played by classical chaos 
was provided. The sensitivity of dynamical tunneling to 
the underlying phase space structure can have important 
ramifications in the molecular context. Phase space of 
molecular systems at higher energies, corresponding to 
high overtone transitions, are typically a mixture of regu- 
lar and chaotic regions. The presence of stochasticity can 
lead to the dynamical tunneling being enhanced (or sup- 
pressed) by several orders of magnitude and hence highly 
mixed states and complicated spectral patterns. Even in 
the near-integrable regimes tiny induced resonances can 
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arise as a result of the primary resonances which can 
enhance or suppress dynamical tunneling between sym- 
metry related modes. Alternatively, molecule-field inter- 
actions can2i lead to creation (destruction) of nonlinear 
resonances which can lead to significant enhancement or 
supression of energy flow. 

The purpose of this work is to establish the role of clas- 
sical resonances in the near-integrable phase space for dy- 
namical tunneling and the nontrivial effect of higher or- 
der induced resonances and perhaps chaos. We focus on 
the symmetric case in this work but the role played by the 
nonlinear resonances is expected to hold in the nonsym- 
metric situations as well. In the near-integrable regime it 
is possible for more than one nonlinear resonance zone to 
manifest in the phase space. In such instances the single 
rotor picture^°!-2i is not sufficient to account for dynam- 
ical tunneling. Yet it is shown that hopping across the 
various resonance islands is an approximate picture that 
yields qualitative (and perhaps quantitative) insights into 
the splitting patterns. We start by motivating the model 
Hamiltonian in section [Til In section IIIII the concept of 
dynamical tunneling and quantum viewpoints are illus- 
trated. The phase space viewpoint and relations to the 
quantum viewpoint are discussed in section HVl bv focus- 
ing on specific cases. In section HVBI two issues regard- 
ing avoided crossings and influence of very small induced 
resonances are addressed briefly. Complications due to 
multiple resonances in the near-integrable limit and the 
resulting resonance-assisted tunneling are illustrated in 
section IIV CI The possibility of chaos-assisted tunneling 
occuring in the system is explored in section IIVDI A 
brief summary is provided in section fVl 



II. MODEL HAMILTONIAN 

As discussed in the introduction attempts to under- 
stand dynamical tunneling, from a classical phase space 
viewpoint, in three or more degrees of freedom systems is 
ambitious without studying the detailed correspondence 
in lower dimensions with multiple resonances. In this 
regard the effective spectroscopic Hamiltonian for water 
due to Baggott-SS provides a good model system. The 
Hamiltonian is given by 

J^C - J^O + g y 1:1 +7^2:2 ^^^f*^ 2:1 +^^2:1 ) l'^) 

with, g' — g + A'(ni + n2 + 1) + X"nf, and 

Ho = ujs{ni + n2) + LUbUb + Xs{nl + nl) 

+ Xbul + Xsbnb{ni + n2) + Xssnin2 (2) 

describing the anharmonic local stretches (1,2), and bend 
(b). The various parameter values are LOg = 3885.57, Wb = 
1651.72, = -81.99,2:6 -18.91, x^f, = -19.12 and, 
Xss — —12.17 cm^^. The j*^ mode occupancy is rij — 
OjUj with (a|,aj) being the harmonic oscillator creation 



and destruction operators. The perturbations, anhar- 
monic resonances, 

K^:p = i4n^,r + i<^lrM' (3) 

connect zeroth-order states |n), |n') with \n'^ — rn] = q 
and \nj — rijl — p. The strengths of the resonances 
are g = -56.48, A' = 3.02, A" = -0.18,7 = -0.91, 
and P — 26.57 cm~^. In order to perform detailed 
studies on the above Hamiltonian in this work we also 
analyze subsystems obtained by retaining specific reso- 
nances. These subsystems, denoted by A and B, are de- 
scribed by the Hamiltonians Ha = Hq + ^{¥2^1^ + ¥2^1^) 
(12) 

and Hb = Ha + 7Vf.2 respectively. The full system 
with all the resonances will be denoted by C. The reason 
for this specific choice of subsystems has to do with the 
fact that classically Ha and Hb exhibit near-integrable 
dynamics whereas He has mixed regular-chaotic dynam- 
ics. H is effectively two dimensional due to the exis- 
tence of the conserved polyad P — {ni + 712) + nb/2 i.e., 
[H, P] = 0. Spectroscopic Hamiltonians can be obtained 
by a fit to the high resolution spectra or from a pertur- 
bative analysis of an existing high quality ab initio po- 
tential energy surface^. In any case such Hamiltonians, 
and variants thereof, have proved quite useful in studying 
the highly excited vibrations of many small molecules'^—. 
Couplings to rotations and large amplitude modes are 
certainly importanti^ but a detailed phase space analy- 
sis of such general rovibrational Hamiltonians is beyond 
the scope of this work. 

The classical limit'^^ of the above Hamiltonian can be 
obtained by using the Heisenberg correspondence, aj <-> 
^yjj exp{i9j), and results in a nonlinear multiresonant 
Hamiltonian 

H{I,e) = Ho{I) (4) 

+ g'^^/hr2Cos{ei - 62) + -fchh cos(26li - 26*2) 
+ f3Jb[^/h cos{ei - 20fc) + /^cos(02 - 26*6)] 

with (Ij,9j) corresponding to the action- angle variables 
associated with the mode j. Note that the actions I 
and the zeroth-order quantum numbers n are related by 
the correspondence I — > (n -t- l/2)?i for simple vibra- 
tions. The various resonant couplings are realted to the 
quantum coupling strengths by g'^ = 2g', 7^ = 27, and 
Pc = 2/?. We note that another reason for the choice of 
this system has to do with the fact that fairly detailed 
classical-quantum correspondence studies of the highly 
excited eigenstates have been already performed'^'^'^^. 

In this study the focus is on a specific set of molecular 
parameters (provided above) which are representative of 
the H2O molecule^i^^. However, at the outset, we em- 
phasize that this choice is by no means special for the 
analysis and conclusions of this work. The goal of this 
work is not to provide accurate quantitative estimates 
for the splittings - something that can be trivially calcu- 
lated quantum mechanically for the model Hamiltonian. 
The emphasis here is on establishing and understanding, 
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FIG. 1: A portion of the energy levels of the model effective 
Hamiltonian corresponding to polyad P = 6. The full system 
is denoted by C and some of the experimental energy levels 
are shown as open circles. B and A show the energy levels 
for the subsystems with specific resonances removed. Some 
of the closely spaced doublets are indicated by *. In C the 
box around the top three levels indicates levels with 4, 5, and 
6 quanta of excitation in the stretch mode. 



qualitatively, the important role of nonlinear resonances 
in dynamical tunneling and for this purpose the above 
model Hamiltonian provides a good paradigm. In par- 
ticular we choose the states corresponding to the polyad 
P = 6 which lie in a energy range of [16600, 21300] cm^^ 
above the ground state. There are two main reasons for 
such a choice. Firstly the motivation is to understand 
aspects of dynamical tunneling for excited vibrational 
states both in the presence of nonlinear resonances and 
chaos^. Secondly the exhaustive compilation'^^ of the ex- 
perimental energy levels of H2O by Tennyson et.al. does 
show three doublets in this energy region. In local mode 
notation these doublets are {60)^0, (50)=^2, and {51)^0. 
A few other local modes are also reported with only one 
of the parity states like (40)~4. Although the effective 
Hamiltonian used here is perhaps not the best it still pro- 
vides some correspondence to the patterns observed in 
the experimental results. Moreover at such high levels of 
excitation, certainly for the bending mode, there can be 
significant interpolyad mixings— which require analyzing 
a three degree of freedom system. Thus we use the ef- 
fective Hamiltonian as a toy model which would guide 
further, more rigorous, studies in full dimensionality. To 
give an idea about the validity of the effective Hamilto- 
nian in the chosen energy range, in Fig. [T] we show the 
energy levels as compared to the experimental levels and 
the agreement is satisfactory. 



III. DYNAMICAL TUNNELING: 
ILLUSTRATION AND QUANTUM VIEWPOINT 

In the rest of the paper the zeroth-order states will be 
denoted by \n1n2) since the bend quantum number rib is 
fixed by the polyad P. In order to illustrate the concept 
of dynamical tunneling we consider zeroth-order degen- 
erate states 1 20) and |02) with only the 2:1 resonances 
present i.e., 5' = 7 = in cq. [TJ As (20|1>2:"^ |02) = = 

(201^2^.^^^ 1 02) the state |20) has no direct coupling to the 
symmetric counterpart |02). Nonetheless there is an indi- 
rect coupling via the 2:1 resonances. However, as shown 
in Fig. [2^ the zeroth-order state is quite far from the re- 
gion of state space which can potentially come under the 
influence of the 2:1 anharmonic resonances. A classical 
trajectory initiated with initial conditions corresponding 
to the state |2,0) will continue to stay in that region of 
the phase space forever without reaching the symmetric 
region corresponding to the state |02). One such classi- 
cal trajectory, for about 0.25 ns, in phase space is shown 
in Fig. [2Jd (inset). This classical trapping occurs despite 
no apparent energetic barriers since in Fig. [5^ we show 
the classical actions satisfying Hq ~ -E20 which connects 
the two states. However, contrary to the classical ob- 
servation, in Fig. [2]3 the quantum survival probability of 
the state |20) shows transfer of population to |02) with 
a period of about 0.20 ns. In a two-state approximation 
this corresponds to a splitting of the degenerate levels 
A2 ~ 0.17 cm~^. This is the phenomenon of dynami- 
cal tunneling. Previous discussions^ ^a^iM of dynamical 
tunneling have used the effective Hamiltonian with only 
the primary 1:1 present. Here we delibrately show this 
effect with the 2:1s only to emphasize the role of induced 
resonances. 

A possible explanation of the tunneling arises from 
the vibrational superexchange^^i^l perspective wherein 
zeroth-order states coupled locally by the 2:1 perturba- 
tions to 1 20) and |02) are considered. One then constructs 
all possible perturbative chains connecting the states |20) 
and |02). An example of such a chain, shown in Fig. [2^, 
is |20) |10) |11) |01) |02). The contribution 
to the splitting is 



f3- 



,4(20|y(i^)|10)(10|t>(2b)|ii)(ii|y(ifc)|0l)(01|y(2'')|02) 



with AE^^^ EE E% 



E° 



(5) 

In principle there are an 
infinite number of chains that connect the two degener- 
ate states. In practice, due to the energy denominators it 
is sufficient to consider the minimal length chains^^. By 
minimal length chains one means inclusion of all possible 
perturbative terms at the lowest order of perturbation 
necessary. For instance in the example above it is neces- 
sary to go to atleast fourth order in perturbation theory, 
thus scaling like /9^, to get a finite value for the splitting. 
In our case there are six minimal chains and summing 
the contributions one obtains a splitting of about 0.18 
cm~^ which compares well with the exact splitting. 
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More generally, splitting for any set of degenerate 
states |rO) and |0r), at minimal order, can be calculated 
as^: 

K^Y.P'^a'l'll^riT^al) (6) 
ahc fi 

with being an index of all possible chains Tabc for 
a particular choice of a, b, c satisfying the constraint 
a + 2b + Ac — 2r. Indeed such a calculation can be 
done on our model system and the resulting perturba- 
tive estimates are compared to the exact splittings in 
Fig. [3] For further discussions the superexchange calcu- 
lations have been compared to the exact results for cases 
A, B, and C. At the outset note that the perturbative 
calculations essentially reproduce, to quantitative accu- 
racy, almost all of the splittings. This degree of accuracy, 
given high orders of perturbation and number of chains, 
is pleasantly surprising. For instance the inset to Fig. [3] 
shows that 130 contributions of varying signs, some as 
large as 10~^, conspire to yield the correct splitting of 
about 10~^ for Aq. Two important observations can be 
made at this juncture. First, the significant role of the 
tiny 2:2 resonance is clear in that addition of this 2:2 to 
the 2:1s enhances the splittings by many orders of mag- 
nitude - something which is not readily apparent from 
Fig. [TJ For even values of r, especially r — 2, this is 
expected since the 2:2 provides a dominant coupling. In 
fact for r = 2 one expects A2 ~ 47 = 3.6 cm~^ which 
compares very well to the exact value of 3.8 cm~^. How- 
ever, it is significant to note that such a calculation for 
r = 6 is too small by a factor of five and that the odd 
r states would not split at all in the presence of a 2:2 
coupling alone. Secondly in all three cases A,B, and C 
the splittings more or less monotonically decrease with 
increasing stretch excitations. The one exception is in 
the full system (case C) for the state with r = 4. Later it 
is shown that this is an instance of suppression of tunnel- 
ing due to an avoided crossing between the doublet and 
a third 'intruder' state. Incidentally in the experimental 
dataiSi pertaining to r = 4 among the doublets, denoted 
in the local mode basis as (40)^4, only the state (40)^4 
is reported^. This, given our model Hamiltonian, could 
be a pure accident! 

It is important to note that the superexchange ap- 
proach invokes the 2:1 resonant terms without any ref- 
erence to the classical phase space. This surprisingly 
good accuracy seems to hold whenever the phase space 
is integrable or near-integrable and is perhaps related 
to the equivalence of the superexchange to the WKB 
approached. If this observation is true in multidimen- 
sional systems then a phase space perspective on dynam- 
ical tunneling should be able to provide some insights. 
Apart from providing a phase space analog to the su- 
perexchange approach it would also be possible to study 
the extent of sensitvity of dynamical tunneling to the 
various classical structures. The rest of the paper is ded- 
icated to uncovering precisely such a phase space picture. 




FIG. 2: Classical and quantum viewpoints on dynamical tun- 
neling for subsystem A. The states are |02) and |20). (a) 
The state space for polyad P = 6 is shown with the regions 
strongly influenced by the 2:1 anharmonic resonances. An 
example of a minimal (solid) superexchange chain connecting 
the states is shown. The energy contour satisfying Ho ~ -E20 
is shown connecting the two states, (b) Quantum survival 
probabilties for [20) and |02) (squares) are shown for 0.25 ns. 
Note the almost coherent transfer of population with a period 
of 0.20 ns. The inset shows a classical trajectory initiated near 
|20) for a time period of about 0.25 ns. The time variations 
of the actions (/i — 1/2) and {I2 — 1/2) are shown as filled 
and unfilled circles respectively. The actions stay trapped for 
the entire period. 




FIG. 3; Comparison of the exact quantum (symbols) dynam- 
ical tunneling splittings A, in cm~^, to the superexchange 
results (lines). Shown are the splitting between the states 
\ni — r,n2 — 0) and |0r). The r = 1,2 cases are not shown 
for the full system since on inclusion of the large 1:1 they are 
not local modes but resonant normal modes with large split- 
tings. The experimental values for the splitting for r = 5, 6 
are shown as filled circles for comparison. The arrow indicates 
the substantial deviation of the superexchange result from the 
exact value for r = 4. The inset shows the contribution from 
130 terms for the superexchange calculation of Ae in case A. 
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IV. DYNAMICAL TUNNELING: PHASE SPACE 
VIEWPOINT 

A. Induced resonances and dynamical barriers 

In what follows, the classical Poincare surface of sec- 
tions are plotted in ((/i - /2)/2, (6'i - 6*2)) = {Ki/2, 20i) 
at constant values of the sectioning angle 4> = 0i+92 — ^9b 
with (/) > 0. In this representation the normal-mode reso- 
nant regions appear around Ii ~ /2. Local modes regions 
correspond to above and below the normal-mode regions 
and large 2:1 bend/local-stretch resonant islands appear 
in the local mode regions. 

Wc begin by asking the question as to what possible 
structure in the classical phase space is mediating the 
dynamical tunneling in Fig. ^Bp. The answer, based on 
earlieri^i2ii^i22ii^ works, should involve a nonlinear res- 
onance zone in the phase space. Indeed the Poincare 
surface of section at _E = £'20 ~ 18641 cm~^ shown in 
Fig. m indicates a resonance zone juxtaposed between the 
zeroth-order states as seen in the phase space. The sur- 
face of section is consistent with the state space view in 
Fig. [5^ in that the primary 2:1 resonances do not ap- 
pear in the phase space. The 2:1 resonance zones, if 
present, would show up in Fig. [2^ around Ki/2 w ±2.6. 
Hence a direct involvement of the 2:1 resonances is ruled 
out but surely the existence of the 2:1s must be giv- 
ing rise to a dynamical barrier, in the form of the ob- 
served resonance zone, mediating the tunneling. In or- 
der to confirm that this nonlinear resonance is mediating 
the dynamical tuneling it is necessary to extract an ex- 
plicit expression for the dynamical barrier and obtain the 
splitting. For this purpose we apply methods of nonlin- 
ear dynamics^i^ on the classical Hamiltonian eq.4 with 
5 = = 7 and (3^0. The analysis below is not new 
and can be found in many of the earlier works includ- 
ing Ozorio de Almeida^S and SibertS^. For completeness 
sake only essential results are explicitly given and tech- 
nical aspects of classical perturbation theory will not be 
discussed. 

To start with a canonical transformation is performed 
{1,9) using the generating function 



F=i9i- 29b) Ji 



- 29h)J2 



(7) 



The choice of the generating function implies 4'k — 9k — 
29b, Jk = h for A: = 1, 2 and 1^3 = 9b,N = 2{h+l2) + Ib. 
The resulting Hamiltonian 

Hi3,xP;N) = Hoi3;N)+ef],iN -2J1-2J2) (8) 
X [\/Ji cos^/^i -I- v^cos V'2] 

is ignorable in the angle "03 implying the conserved quan- 
tity N = P + 5/2 which is the classical analog of the 
polyad number P. The above two-dimensional Hamil- 
tonian contains the 2:1 resonances and hence noninte- 
grable. However the phase space from Fig. [3] suggests 
that at E' « E'j 12 the 2:1s do not have a direct role 



to play. Another way of stating this is that the nonlin- 
ear frequencies ilk = dHQ{T)/dIk are far away from the 
condition ^1^2 — 2flb = for the states in consideration. 
For the system parameters in this study the |2, 0, 8) state 
is detuned by almost 800 cm^^. Consequently a formal 
parameter e has been introduced with the aim of pertur- 
batively removing the 2:1 resonances, characterized by 
-01,2 in eq. [HI to 0(e). This can be done by invoking the 
generating function 



G = JitI)i + J2'02 + sin -01 -t- 52 sin 02 



(9) 



where the functions 51,2 — gi,2iJi, J2) are determined 
by the condition of the removal of the primary 2:1s 
to 0{e). The angles conjugate to J are denoted by 
i/j. The new variables {3, if)) are related to the vari- 
ables {3,xp) through the relations Ji^2 — 9G/90i_2 and 
01,2 = dG/dJi^2- Using the identities: 



cos(e2; sin((/7)) 
sin(e2; sin(<y9)) 



Jo{ez) + 2^2^21 (ez) cos(2Z^)(10) 
2^^2/+i(ez)sin((2Z+l)</7) (11) 



with J7fe being the Bessel functions, it can be shown that 
H{3,'4,;N) = Ho{3;N) + Y,^''Hki3,'4>;N) (12) 

k 

where 



Ho = ns{Ji + J2) + assiJ! + Ji) + 



with ils — LOs — 2u!b -\- {xsb — 'ixb)N, as 
and ai2 — 8xb — 'ixsb + Xss- At 0{e) 



0112J1J2 

= Xs+^Xb- 



(13) 
- 2xsb 



Hii3,'iP;N)= ^ (/?,JbV^+<?fe(J)f^fe)cos0fe (14) 



fc=l,2 



where we have denoted Jb = N — 2(Ji -I- J2), and (ik = 
dHo/dJk- The action Jb is essentially the number of 
bend quanta in the system. The choice of the fmictions 



5fc(J) = 



Jb 



(15) 



eliminates the 2:1s to 0{e) and hence the O(e^) term in 
the transformed Hamiltonian eq. [T^lis determined to be: 

H2i3,tp;N)= Afe(H-cos20fc)-HA±cos(0i±02) 



fe=l,2 



(16) 

The term cos(0i — 02) is easily identified with a 1:1 res- 
onance between the modes 1 and 2 and the strength of 
this induced 1:1 resonance is obtained as 



A_(J;iV) = i/32 



2{ni + n2) + ai2Jb 
nih2 



JbVjJ'2 (17) 
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FIG. 4: Classical Poincare surface of section for Ha at _E ~ 
i?20 showing the induced 1:1 resonance island due to the two 
stretch-bend 2:1 resonances. 



Expressions for the strengths Ai_2 and A+ are not needed 
for our analysis hereafter and hence not explicitly given 
in this work. 

At this stage the transformed Hamiltonian to O(e^) 
in eq. [H] still depends on both the angles V'l and ip2 
and hence non-integrable. In order to isolate the in- 
duced 1:1 resonance we perform a canonical transforma- 
tion to the variables (K, 0) using the generating function 
G = (Vii - 'ip2)Ki/2 + (V^i ip2)K2/2 and average the 
resulting Hamiltonian over the fast angle 02- The reso- 
nance center, ~ 0, approximation is invoked resulting 
in a pendulum Hamiltonian describing the induced 1:1 
resonance island structure seen in the surface of section 
shown in Fig. Within the averaged approximation the 
action K2 = Ji -I- J2 is a constant of the motion and can 
be identified as the 1:1 polyad associated with the sec- 
ondary resonance. The resulting integrable Hamiltonian 
is given by 



H{Ki,^uK2,N) 
where 

g,nd{K2,N) 

with 



1 



2Mi 



-Kf + 2g,nd{K2,N) cos 201 



2(q!i2 - 2assy 

0, 



(18) 



(19a) 



^f{K2,N){N ~2K2)K2 (19b) 



KK2,N) 



4(»3 + assK2) + ai2N 
[2(ns + assK2) + ai2K2]' 



(19c) 



Note, and this is important for the discussion in sec- 
tion IIVBI the averaging procedure does not remove the 



FIG. 5: Variation of the |20)± doublet energies with the 2:1 
resonance strength j3. Results correspond to Ha and no exact 
or avoided crossings are observed. The inset compares the 
splittings A2 calculated semiclassically (open circles, eq. I20|) 
to the exact quantum (solid line) over the same range of /3 
variation. Note that A2 increases by a factor of 10** which 
is captured well by the induced 1:1 resonance. All quantities 



are m cm 



higher harmonics of the induced 1:1. Thus by focusing 
on the induced 1:1 alone in eq. [TH] we have neglected 
all induced resonances of the form q : q. In terms of 
the zeroth-order quantum numbers Ki = ni — n2 and 
K2 = ni -I- 712 + 1 = m + 1. Note that eq. I19bl more 
generally eq. 1171 is an analytic expression for the dynam- 
ical barrier seperating any two localized states in terms 
of the zeroth-order actions (/1,/2,/b). 

One can now use eq. [15] to calculate the dynamical 
tunneling splitting of the degenerate modes \ni = r, ^2 = 
0, Ub = 2{P-r)) and \ni = 0, 712 = r, 7ib = 2(P-r)) via^i 



9i: 



(r-2) 

' n 



9i7id 



El{r)-El{k) 



^™'^2'^-i[(r-l)!]2 



(20) 



where E^{k) — fc^/2A/ii is the zeroth-order energy. For 
our example with r = 2, 771 = 2 using the parameters 
of the Hamiltonian we find Mn w 1.32 x 10~^ and 
gind ~ 3.93 cm^^. Within the pendulum approximation 
the half-width of the resonance zone can be estimated 
as 2^/2Miigind ~ 0.64. This estimate, relalizing that 
Ki « (/i — 12), is in good agreement with the surface of 
section in Fig. [4] The resulting splitting Aj'^ « 0.20 cm~^ 
agrees well with the exact splitting. This proves that the 
induced 1:1 resonance arising from the interaction of the 
two primary 2:1 resonances is mediating the dynamical 
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tunneling between the degenerate states. Moreover, from 
a superexchange perspective it is illuminating to note 
that the splitting can be calculated trivially by recogniz- 
ing the secondary phase space structure (a viewpoint em- 
phasized in Ref. 27 as well) . This observation emphasizes 
the superior nature of a phase space viewpoint on dynam- 
ical tunneling. The analysis above also sheds light on the 
relation between the quantum superexchange viewpoint 
and WKB approximation. Note that gind ^ 0^ and thus 
the phase space viewpoint, with A2 ~ gf^^^ is equivalent 
to the superexchange approach with A2 ~ In this 
near-integrable situation this correspondence generalizes 
the observation by Stuchebrukhov and Marcus about the 
equivalence of superexchange and WKB methods for in- 
tegrable systems^. 

B. Higher harmonics and level crossings 

In this section we address two issues relevant in the 
context of dynamical tunneling. The first issue has 
to do with the role and relation of eigenvalue avoided 
crossing o^^'^^ to nonlinear resonances and hence the phe- 
nomenon of dynamical tunneling. Since we have just ar- 
gued for the involvement of the 2:1 resonances in medi- 
ating the dynamical tunneling between degenerate states 
it is interesting to ask if the states are also involved in 
an avoided crossing as the 2:1 strength /3 is varied. The 
answer to this question is negative and Fig. [5] provides 
evidence that the states |20)± do not undergo any iden- 
tifiable avoided or exact crossings with varying f}. In fact 
the splittings (Fig. [5] inset) predicted from eq. [20] agree 
fairly well over the entire range of /3. It is also worthwhile 
noting that for the system Ha the splittings Ai.2,3 are 
described quite well by a single induced 1:1 resonance. 
Therefore, in general, it is not necessary that dynami- 
cal tunneling between two zeroth-order states \n1n2) and 
|ni n2 — q) be mediated by a p : 5 nonlinear resonance 
alone. 

A second issue has to do with the role of higher har- 
monics or Fourier components in dynamical tunneling^^ . 
The doublet splitting A2 for Ha was seen to arise from 
a induced 1:1 resonance. The perturbative analysis, on 
the other hand, even in the averaged approximation con- 
tains all the higher harmonics of the induced 1:1. These 
Fourier components appear at higher orders in e and have 
very small, but finite, contribution to A2. If for some rea- 
son the leading order induced 1:1 were to vanish then the 
higher Fourier components become important. For the 
present example the 2:2 component would be especially 
significant. In order to demonstrate this effect consider 
the Hamiltonian 

H^Ho- \g\vl^^^ -t- P{V^^^ + V^^^) (21) 

A 1:1 resonance has been added to Ha with a strength 
opposite in sign to the induced 1:1 resonance. The 
strength g is varied for fixed /3 and there is a specific 
value of 5 = go for which the primary and the induced 



1:1 resonance cancel each other out. The value of g^ can 
be estimated by recognizing the integrable pendulum ap- 
proximation to the above Hamiltonian as 

H^T;^Kl + [-\g\K2 + 2g,^d)cos2(^i (22) 

Thus around jf/jo ~ 2gind/K2 the above leading order 
approximation would predict a shutdown of dynamical 
tunneling. From the estimates provided earlier it is easy 
to obtain the value |g|o ~ 2.62 cm~^ for A2. 

In Fig. [6^ the variation of the exact splitting A2 is 
shown as a function of g. The results clearly show the 
destruction of tunneling near go but a more significant 
observation is that at goj labeled by region ' M' in Fig.[6K, 
the actual value of the splitting is far from zero. This is 
in contrast to the prediction of shutdown of dynamical 
tunneling based on the integrable pendulum approxima- 
tion above. The corresponding phase space in fig. [6lVI 
indeed shows that the 1:1 resonance zone has almost van- 
ished. So what is mediating the dynamical tunneling in 
region Ml Arguments based on some residual 1:1 cou- 
pling are quickly dispelled by noting that the width of 
the resonance zone in Fig. [HM is nowhere close to the 
required value. Chances of a strong avoided crossing 
leading to the enhanced splitting are ruled out by in- 
specting Fig. Perhaps there is a broad avoided cross- 
ing, as suggested by the exact crossing of the doublets 
|20)± in fig. (inset), leading to the two dips in A2 
(cf . fig. [6^) , with the upper state which happens to be 
a normal mode state. It is relevant to note that the en- 
ergy seperation between the doublet and normal mode 
state is almost equal to the mean level spacing. More- 
over a superexchange calculation of A2 at M including 
the required minimal chains (cf . eq. [5]) r4oo , ro2o , and 
r2i0j yields a value 1.70 x 10"'^ cm~^ which is in excel- 
lent agreement with the exact value of 1.65 x 10"'^ cm~^. 
All this suggests that in the region M a tiny, but leading 
order, 2:2 with strength 7 — A2/4 « 4.1 x 10"'* cm"^ is 
mediating the dynamical tunneling. A few important ob- 
servations provide further support and we briefly outline 

them. Firstly the splitting in region M are very sensitive 

(12) 

to adding a perturbation ±71/2.2 to the Hamiltonian in 
eq. [21] whereas the regions L and R are far less sensitive 
as shown in Fig. [S^. Secondly, the classical perturbation 
analysis carried out to higher orders reveals the emer- 
gence of an induced 2:2, scaling as (i'^g at O(e^), with 
strength comparable to 7. Indeed the exact A2 in region 
M for different values of (3 scales linearly with pP'g. The 
superexchange calculation also shows that the contribu- 
tion from the minimal family r2io *.e., 0^g nearly, but 
not exactly, balances the contribution from and g^ 
families. The analysis here is an example of what Pear- 
man and Gruebele^ call as 'phase cancellation' effect in 
perturbative chains. Evidently, the analog of this effect 
in phase space has to do with the destruction of nonlin- 
ear resonances. Although only the r = 2 case has been 
shown here, similar effects arise for the higher excitations 
and can be more complicated to interpret. 
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FIG. 6: (a) Variation of the exact A2 with the primary 1:1 
strength g. Note the two dips and the relatively flat region 
denoted by M. The arrow indicates the strength go. The 
effect of adding a small 2:2 with strength 7 ~ ±4.1 x 10"'* 
cm~* at specific places are shown as triangles. Inset shows 
the linear scaling of exact A2 (points) with /S^I^Iq. In (b) the 
variation of the relevant energy levels are shown. Inset shows 
the details for the doublet |20)± which undergo two exact 
crossings. The three panels denoted by L, M, and R show 
the surface of sections corresponding to the regions labeled 
in (a). The sections are computed using the full classical 
limit Hamiltonian corresponding to Ha — \g\Vi^-^^ ■ Note the 
reorganization of the phase space in going from L to R. 



C. Multiple resonances: Resonance-assisted 
tunneling 

One of the key lessons learnt from the preceeding dis- 
cussion is that nonlinear resonances have to be identified 
in order to obtain a clear phase space picture of dynam- 
ical tunneling. Given the near-integrable phase space of 
Ha it is tempting to apply eq. [201 to calculate the split- 
ting of the other doublets corresponding to increasing 
stretch excitations. In particular the monotonic decrease 
with increasing stretch excitations, seen in Fig. [31 seems 
to fit well with the notion of a single resonance (the in- 
duced 1:1 here) mediating the tunneling. The naivety 
of such a viewpoint is illustrated in Fig. [7^. Clearly for 
r = 1,2,3 the expectation holds but deviations arise al- 
ready at r = 4 and the splittings are many orders of 
magnitude smaller than the quantum results for r — 5,6. 
The reason for this large error is immediately clear from 
the surface of sections shown in Fig.[7lD,c,d. The primary 
2:1 resonance has appeared in the phase space and coex- 
ists with the induced 1:1. The states corresponding to 
r = 4, 5 are 'resonant' local modes and a Husimi repre- 
sentation of the eigenstate in the classical phase space (cf. 
Fig. [TUk ) confirms their nature. The r = 5 doublets are 




2 4 6 




0.5 0.75 1 

FIG. 7: (a) Comparison of the exact (open triangles), su- 
perexchange (solid line) , and the phase space resonance based 
(filled triangles) splittings for Ha versus excitation quanta in 
the OH-stretch. Multiple resonance corrections to r = 5, 6 
are shown by shaded circles. See text for details. In (b), (c), 
and (d) the surface of sections at energies corresponding to 
r — 4,5, and 6 respectively are shown. The axes ranges for 
(b) and (d) are identical to those shown in (c). In (d) one 
possible dynamical tunneling sequence is shown by arrows. 



in fact localized in the large 2:1 resonance islands. This 
can also be anticipated from the state space location of 
the resonance zones shown in Fig. [5^. Thus a single ro- 
tor integrable approximation is insufficient to describe 
the dynamical tunneling for such cases. Notice the rich 
structure of the phase spaces in Fig.[7)D,c,d resulting from 
the presence of two 'distant' and 'nonoverlapping' (in the 
Chirikov^ sense) resonances - a clear manifestation of 
the nontrivial nature of the near-integrable systems. 

Is it possible to use the phase space structure, for ex- 
ample in Fig. [Tli, and calculate the splitting Ag? The 
phase space shows the existence of the primary 2:1s, the 
induced 1:1, and a thin multimode resonance. For now we 
will ignore the multimode resonance and consider the pri- 
mary 2:1s and the induced 1:1 to be independent. This is 
a strange approximation given that both the induced 1:1 
and the multimode arise due to the 2:1s. However, since 
they are 'isolated' from each other and lacking any other 
simple approach we take this viewpoint and compute 
the splitting. This approach has been advocated ear- 
lier by Brodier, Schlagheck and UUmo in their analysis— 
of resonance- assisted tunneling in kicked maps. In the 
present case the effective h is too large to be in the semi- 
classical limit and hence quantitative accuracies are not 
expected. At the same time this is a blessing in dis- 
guise since we do not have to deal with the the myriad 
of resonance zones in the phase space! Schematically we 
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imagine the following (cf. Fig.[7li): 



|60) ^ |40) ^ |04) 



|06) 



(23) 



where (3eff is the effective coupling across the 2:1 islands 
and gind is the effective 1:1 coupling derived in the pre- 
vious section. It is important to note that the gind value 
appropriate to r = 6 needs to be used i.e., the resonance 
zone width as seen in the surface of section in Fig. [6ji. 
The effective 2:1 couplings can be extracted by approxi- 
mating the 2:1 resonances, for instance the top island in 
Fig. [Tji, by a pendulum Hamiltonian: 



1 



2M2 



-(Ji - Jl)^ + 2l3pCosiPi 



(24) 



withMzi = {2ass)~^ and/3p = /?c(iV-2 J2-2J[)yjf/2. 
The resonance center J[ = —M2i{2Qs + Q!i2)/2 « 5.7. 
The effective coupling is thus estimated as: 



^6 



El 



(25) 



with E° = {h{ni + 1/2) - J[)V(2Af2i). The splitting is 
then calculated via the expression 



A. 



El -El 



(26) 



where Al'^ is the splitting for r — 4 calculated assuming 
the induced 1:1 resonance alone. Using the numerical val- 
ues for the various parameters we find Ag ~ 4.8 x 10^^ 
cm^^ as compared to the exact value of 1.4 x 10~^ cm~^. 
This simple estimate is quite encouraging if one also con- 
siders the fact that the induced 1:1 alone would give 
~ 10~^^ cm~^ and the superexchange value is 6.5 x 10"'' 
cm~^. Similar calculation for r = 5 doublet also im- 
proves the result as shown in Fig. [7^. This simple minded 
scheme has been tested for other parameters and the re- 
sults are comparable to the present case. Further support 
for such an approach comes from considering the r — 6 
doublet splitting with an additonal 2:2 resonance present 
i.e., described by the Hamiltonian Hb. From Fig. [3] it 
is clear that the addition of the 2:2 results in about four 
orders of magnitude increase in Ag. The phase space is 
still near-integrable and a calculation based on the 2:2 
alone underestimates the exact result of 3.7 x 10~^ cm~^ 
by nearly a factor of four. Once again employing the 
approach of hopping across the 2:1 using (3eff, followed 
by the 2:2 connecting |40) to |04), and hopping across 
the symmetric 2:1 island improves the result yielding a 
splitting of about 2.8 x 10^"^ cm"-'^. Interestingly the su- 
perexchange calculation has a dominant contribution of 
similar magnitude from the family r402 which is readily 
identified with the /3'*7^ 'path' in phase space outhned 
above. The analysis here indicates that it is conceivable 
that nonlinear resonances can lead to dynamical tunnel- 
ing between energetically similar regions of phase space 
even in nonsymmetric situations. 
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FIG. 8: Variation of eigenenergies with the 1:1 strength pa- 
rameter g. (a) r = 4 doublet involved in an avoided crossing 
with a — parity normal mode state. Inset shows the behaviour 
of A4 with varying g. (b) Same as in (a) for the r — 5 dou- 
blets. All parameters are in cm~^. Note that the inset x-axis 
range for g is from to —60 cm~^. Chaos is visible in the 
phase space for \g\ > 30 cm^^. 



D. Mixed phase space: Chaos-assisted tunneling? 

Up until this point dynamical tunneling was investi- 
gated with the underlying phase space being integrable 
or ncar-integrable. We now look at the full system 
He which exhibits mixed chaotic-regular phase space. 
The phenomenon of chaos-assisted (suppressed) tunnel- 
ing (CAT) has to do with the coupling of quantum 
states, localized on two symmetry-related regular regions 
of phase space, with one or more 'irregular' states delo- 
calized over the chaotic sea3^. Clearly such processes 
cannot be understood within a two-state picture and the 
sensitivity of the nature of the chaotic states to paramet- 
ric variations is manifested in the splitting fluctutations 
of the tunneling doublets. It is also reasonable to expect 
that the effective h in the system needs to be sufficiently 
small in order for CAT to manifest itself. Perhaps it is 
useful to note a recent analysis^ by Mouchet and De- 
lande wherein they argue that even in cold atom tunnel- 
ing experiments the effective h is not small enough to 
observe CAT. The system investigated herein is nowhere 
close to such limits and at the outset we do not ex- 
pect to see CAT. In other words for our full system de- 
scribed by He only resonance-assisted tunneling (RAT) 
mechanisms should be sufficient to explain the relevant 
splittings^. 

We focus on the doublets r = 4, 5, 6 since the phase 
space shows substantial chaos at the corresponding en- 
ergies. In Fig. [3] it is seen that the r = 4 case showed 
marked deviations from monotonicity. Moreover the su- 
perexchange calculation predicted a much higher A4. In 
Fig. [5^ we show the variation of energy levels with the 1:1 
strength g and a clear avoided crossing between one of 
the tunnel doublets (with — parity) and a normal mode 
state is observed. The inset in Fig. [5^ implicates the 
avoided crossing in the observed exact crossing between 
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FIG. 9: Husimi distributions of the r — 4 doublets and a 
close lying normal mode state for the full system He- The 
axes ranges are identical to those shown in fig. [TId. The dark 
regions correspond to large values while the light regions sig- 
nify smaller values. The top row shows the Husimis for the — 
and + parity states while the bottom left is for the -- parity 
normal mode state. The bottom right panel shows the phase 
space with clearly identified 1:1 and 2:1 resonance islands. 
Note the delocalized Husimis of the doublets in contrast to 
that of the normal mode state. 



the tunnel doublets themselves leading to the suppression 
of tunneling. In F ig, p the phase space Husimi^^ repre- 
sentations (refer to [33 for details regarding the computa- 
tion of the Husimi distributions) of the three states are 
shown with the corresponding phase space. Notice that 
the normal mode ('intruder') state has a much more regu- 
lar Husimi in the phase space as compared to the doublet 
Husimis themselves. In particular the — doublet is exten- 
sively delocalized over the phase space. Given the value 
of the effective h it is not possible to declare the — par- 
ity doublet to be a chaotic state. Note however that the 
doublets corresponding to r = 4 are essentially living in 
the chaotic regions. For the near-integrable subsystems 
Ha and Hb the r = A doublets live on the separatrix 
associated with the 2:1 resonance. So this is an instance 
where a regular intruder state is affecting the tunneling 
between irregular doublets. No attempt is made here to 
calculate the tunneling splitting semiclassically. However 
we note that the exact value A4 « 0.1 cm~^ as compared 
to the superexchange value of 2.3 cm~^. Interestingly 
with further increase of g, resulting in more chaos, the 
superexchange calculation is only a factor of two higher 
than the exact value. 

On the other hand in Fig. [81d the variation of the r = 5 
doublet energies are shown as a function of the 1:1 pa- 
rameter g. This case does not show any avoided cross- 
ing. The Husimis for this doublet corresponding to the 



systems Ha,Hb (near-integrable) and He (mixed) are 
shown in Fig. [TUk.b. and c respectively. It is notable 
that the splitting in going from Ha to Hb increases by 
more than three orders of magnitude (cf . Fig. ^ . The 
Husimis and surface of sections in Fig. [TOk.b are quite 
similar except for the 2:2 resonance zone clearly visi- 
ble in Fig. [TUb . However, the 2:2 alone cannot mediate 
dynamical tunneling between the localized states in the 
2:1 zone. As in the previous section, excitation out to 
the seperatrix state |40) followed by tunneling across the 
2:2 zone to |04) and finally coupling to the symmetric 
1 05) state captures quantitatively the tunneling enhance- 
ment. This is a clear manifestation of RAT in the sys- 
tem. Adding the primary 1:1 resonance, resulting in He, 
the doublet splitting is further enhanced by two orders 
of magnitude (cf.Fig. The phase space as shown in 
Fig. [To] is clearly of a mixed nature. The large 1:1 res- 
onance zones are visible and tentatively one assigns the 
enhancement in splitting to the 1:1 coupling. The 1:1 
alone can account for one order of magnitude enhance- 
ment. However of all the possible resonant subsystems 
considered, quantum mechanically or semiclassically, the 
best estimate is a factor of two too small as compared 
to the A5 for the full system. In this context an inspec- 
tion of Fig. [TOj Husimi for the full case reveals amplitude 
spreading through the classically chaotic region. This 
is in stark contrast to the well localized Husimis in the 
near-integrable cases. Again the effective Ti is too large 
to implicate the classical stochasticity for part of the en- 
hancement in A5. Although the case for r = 6 is not 
shown here Ag is enhanced (cf. fig. [3]) by four orders of 
magnitude in going from Ha to Hb- This enhancement 
could again be captured quantitatively by the semiclassi- 
cal analysis. However the best integrable estimate based 
on the 1:1 and the 2:2 resonances is almost an order of 
magnitude too small for the full system He- 



The upshot of the precceding analysis is that no in- 
tegrable subsystem can account for the splittings of He 
and one atleast needs some chaos to be present. Even 
subsystems with chaos, for instance 2:1s -I- 1:1, on exact 
diagonalization yield spittings which are about a factor 
of two in error. This points to a fairly subtle involvement 
by all the resonances, irrespective of their strengths, and 
is not understood well at this point of time. Reducing 
the effective h in the system is a theoretical tool to clar- 
ify the picture. For our system, a nonscaling one, it is 
not very easy to perform ?i-scaling computations. Nev- 
ertheless some preliminary calculations indicate that the 
near-integrable cases show exponential fall off with h'^^ 
whereas the mixed cases exhibit fluctuations. In par- 
ticular the r = 4 doublet splitting exhibits an algebraic 
dependence h~°' rather than the usual exp(— l/?i) depen- 
dence. 



12 




FIG. 10: The top panels show the Husimis superposed on 
phase space for the + symmetric partner of the r = 5 doublet 
in the near integrable cases. The axes ranges are identical to 
those shown in fig. [7J;. Top left is for Ha (2:1 only) and top 
right is for Hb = Ha + 2:2. The Husimis in both instances 
are localized in the 2:1 islands. The bottom left panel shows 
the Husimi for the full case He and the corresponding phase 
space is shown in the bottom right. Note the small isthumus 
across the chaotic sea connecting the two islands. 



V. CONCLUSIONS 

In this work an intimate connection between dynami- 
cal tunneling and the resonance structure of the classical 
phase space is established. The order and widths of the 
resonances determine the explicit dynamical barriers in 
the zeroth-order action (quantum number) space. The 
current study indicates that it is safer to 'blame' dy- 
namical tunneling on the nonlinear resonances without 
invoking further connections between avoided crossings 
and the resonances. If accidental avoided crossings oc- 
cur, as they do in near-integrable, mixed and chaotic 
systems, then dynamical tunneling can be enhanced or 
suppressed further. This study also lends support to ear- 
lier suggestions^^ that even very small resonant couplings 
can be the cause of the experimentally observed narrow 
spectral clusters. The resonances responsible for such 
spectral clusters can be primary or induced resonances 
and with small but significant amount of stochasticity in 
the phase space the consequences for energy flow can be 
nontrivial. 

In near-integrable situations, generic to lower energy 
regimes in large molecules, a combination of more than 
one resonance zone can control dynamical tunneling. Not 
only is it important to identify the relevant resonance 
zones but it is also necessary to prescribe a route to cal- 
culate the splittings. It is significant to note that nei- 
ther of these problems are easy to solve. Again, and 



not surprisingly, the lack of understanding of the struc- 
ture of multidimensional phase space is the main bottle- 
neck. In this regard the superexchange approach works 
well but such an accuracy for high dimensional systems, 
with mixed regular-chaotic phase spaces, is not guaran- 
teed. As shown in the present work and in a previous 
work'^^ the superexchange approach is prone to errors, 
atleast by an order of magnitude, whenever states are 
involved in avoided crossings. It remains to be seen if 
going beyond the minimal paths prescription, with some 
clever 'resummation' tricks, would lead to quantitative 
improvements. We have shown that treating the res- 
onance zones as pendulums and coupling across them 
to connect the two nearly degenerate states provides an 
approximate route to calculating the tunnel splittings. 
This chain of resonance zones, in some sense, provides a 
semiclassical basis for the success of the vibrational su- 
perexchange approach. The current work has shown the 
close correspondence between various terms contributing 
in the superexchange approach and the underlying res- 
onances in the phase space. The superexchange mech- 
anism itself has its roots in the field of electron trans- 
fer in molecular systems^. In this context it is inter- 
esting to note the analogy between the role played by 
the 'bridge' states connecting the donor and acceptor 
sites for electron transfer— and the resonance zones con- 
necting two symmetry related, distant states localized in 
phase space. Some similarities, perhaps formal, to this 
work showing the promotion of dynamical tunneling by 
multiple resonances and a recent work^^ demonstrating 
promotion of deep tunneling through molecular barriers 
by electron-nuclear coupling further exemplify the close 
parallels between the phenomena of off-resonant electron 
transfer and dynamical tunneling. Regarding the vari- 
ous mechanisms of dynamical tunneling in molecules this 
work suggests that for small molecules with low density of 
states (large effective fi) the dominant mechanism would 
be resonance-assisted. For molecules with sufficiently 
high density of states (sufficiently small effective Ti) both 
chaos and resonances can mediate dynamical tunneling 
and the dominant mechanism is decided by the energy 
of the doublets and the location of the doublets in the 
corresponding phase space. 

Although the analysis was on coupled three-mode sys- 
tems the conclusion remains valid in general for systems 
which exhibit local mode behaviour. This is based on 
our analysis, along similar lines, for many other systems 
which are described more naturally by local (D2O, H2S) 
or normal mode (SO2) limits. For systems with small or 
vanishing anharmonicities, where the issue of the exis- 
tence of localized modes itself is moot, the various reso- 
nances are still important but any calculation based on 
a pendulum approximation is clearly invalid. Thus spec- 
tra of molecules involving light atom stretch-bend modes 
are good candidates to observe the fingerprints of dy- 
namical tunneling. In particular we suggest the high 
overtone spectral regions of the water molecule as one 
possibility. This suggestion is tentative since a model 
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Hamiltonian has been utilized for the analysis and it is 
probable that the polyad picture will break down at such 
high energies. One possibility is to do a detailed analy- 
sis on the best potential energy surface available and/or 
explicitly break the polyad by adding further weak res- 
onances. Both approaches result in genuine three de- 
gree of freedom phase spaces and thus a straightforward 
extension is not easy. Apart from the limitation in vi- 
sualizing the global phase space structures this has to 
do with the fact that there are classical phenomenEsii 
in three or more degrees of freedom that could modify 
certain aspects of the lower dimensional analysis. For in- 
stance in three or higher degree of freedom systems two 
long time phenomena are controlled by nonlinear reso- 
nances - Arnol'd diffusion^^'^^ and dynamical tunneling. 
Precious little is known about either of these phenomena 
and their possible spectral manifestations in such situ- 
ations and needs further study^°. As a final note we 
comment on the possibility of observing long lived local 
excitations in systems with multiple resonances. This 



work shows that additional resonances can modify the 
near-degeneracies via CAT and RAT. But this modifi- 
cation, by the same additional resonances and resulting 
chaos, can go either way due to multistate interactions! 
Moreover, as shown in this work, additional resonances 
and presumably coupling to rotation o^^'^^ can also con- 
spire to increase the degeneracy. Coherent manipulation 
of tunneling using external fields has been demonstrated 
before^ii^ and there are reasons to think that such ideas 
can be utilized in the present context as well. 



VI. ACKNOWLEDGEMENTS 

It is a pleasure to thank Peter Schlagheck for criti- 
cal and illuminating discussions. I am grateful to Prof. 
Klaus Richtcr for the hospitality and support at the Uni- 
versitat Regensburg where part of this work was done. 



^ M. Quack, Annu. Rev. Phys. Chem. 41, 839 (1990); T. 
Uzer, Phys. Rep. 199, 73 (1991); K. K. Lehmann et.al, 
Annu. Rev. Phys. Chem. 45, 241 (1994); D. J. Nesbitt 
and R. W. Field, J. Phys. Chem. 100, 12735 (1996); M. 
Gruebele and R. Bigwood, Int. Rev. Phys. Chem. 17, 91 
(1998); J. C. Keske and B. H. Pate, Annu. Rev. Phys. 
Chem. 51, 323 (2000). 

2 M. E. Kellman, Adv. Chem. Phys. 101, 590 (1997); G. S. 
Ezra, Adv. Class. Traj. Meth. 3, 35 (1998); M. V. Kuzmin 
and A. A. Stuchebmkhov, in Laser Spectroscopy of Highly 
Vtbrationally Excited Molecules, Ed. V. S. Letokhov, p. 
178 Adam Hilger, Bristol, 1989; A. A. Ovchinnikov, N. 
S. Erikhman, and K. A. Pronin, Vibrational- Rotational 
Excitations in Nonlinear Molecular Systems, Kluwer Aca- 
demic/Plenum, New York, 2001. 

3 M. J. Davis and E. J. Heller, J. Chem. Phys. 75, 246 
(1981);E. J. Heller and M. J. Davis, J. Phys. Chem. 85, 307 

(1981) . The notion of associating a generalized tunneling 
with any classically forbidden dynamical event is already 
found in W. H. Miller, Adv. Chem. Phys. 25, 69 (1974). 

* For recent reviews on the local modes see, B. R. Henry 

and H. G. Kjaergaard, Can. J. Chem. 80, 1635 (2002); 

P. Jensen, Mol. Phys. 98, 1253 (2000); L. Halonen, Adv. 

Chem. Phys. 104, 41 (1998). 
^ C. Jaffe and P. Brumer, J. Chem. Phys. 73, 5646 (1980). 
" M. E. Kellman, J. Chem. Phys. 83, 3843 (1985). 

M. Quack, Faraday Discuss. Chem. Soc. 71, 359 (1981). 
® S. Flach and C. R. Willis, Phys. Rep. 295, 181 (1998); J. 

Dorignac and S. Flach, Phys. Rev.B 65, 214305 (2002); S. 

Flach and V. Fleurov, J. Phys.:Condens. Matter 9, 7039 

(1997); V. Fleurov, Chaos 13, 676 (2003) and references 

therein. 

^ R. T. Lawton and M. S. Child, Mol. Phys. 37, 1799 (1979); 
R. T. Lawton and M. S. Child, Chem. Phys. Lett. 87, 217 

(1982) . 

^° E. R. Th. Kerstel et. al, J. Phys. Chem. 95, 8282 (1991). 
A. A. Stuchebrukhov and R. A. Marcus, J. Chem. Phys. 
98, 6044 (1993). 



A. Mcllroy and D. J. Nesbitt, J. Chem. Phys. 91, 104 
(1989). 

A. Callegari et. al, Mol. Phys. 101, 551 (2002). 

A. Portnov, S. Rosenwaks, and I. Bar, J. Chem. Phys. 121, 

5860 (2004). 

T. K. Minton, H. L. Kim, S. A. Reid, and J. D. McDonald, 

J. Chem. Phys. 89, 6550 (1988). 

M. E. Kellman, J. Chem. Phys. 76, 4528 (1982). 

W. G. Harter and C. W. Patterson, J. Chem. Phys. 80, 

4241 (1984). 

K. K. Lehmann, J. Chem. Phys. 95, 2361 (1991); K. K. 
Lehmann, J. Chem. Phys. 96, 7402 (1992); See also, M. S. 
Child and Q. Zhu, Chem. Phys. Lett. 184, 41 (1991). 
R. M. Lees, Phys. Rev. Lett. 75, 3645 (1995); J. Ortigoso, 
Phys. Rev. A 54, R2521 (1996). 
^° A. M. Ozorio de Almeida, J. Phys. Chem. 88, 6139 (1984). 
E. L. Sibert HI, W. P. Reinhardt, and J. T. Hynes, J. 
Chem. Phys. 77, 3583, 3595 (1982. 

J. S. Hutchinson, E. L. Sibert III, and J. T. Hynes, J. 
Chem. Phys. 81, 1314 (1984). 
E. L. Sibert HI, J. Chem. Phys. 83, 5092 (1985). 
'^'^ A. A. Stuchebrukhov and R. A. Marcus, J. Chem. Phys. 98, 
8443 (1993). 

25 E. J. Heller, J. Phys. Chem. 99, 2625 (1995); E. J. Heller, 
J. Phys. Chem. A 103, 10433 (1999). 

R. Roncaglia, L. Bond, F. M. Izrailev, B. J. West, and P. 
Grigolini, Phys. Rev. Lett. 73, 802 (1994); L. Bond, A. 
Farusi, P. Grigolini, and R. Roncaglia, Phys. Rev. E 58, 
5689 (1998). 

O. Brodier, P. Schlagheck, and D. UUmo, Phys. Rev. Lett. 
87, 064101 (2001); O. Brodier, P. Schlagheck, and D. 
UUmo, Ann. Phys. 300, 88 (2002). 
2* O. Bohigas, S. Tomsovic, and D. UUmo, Phys. Rep. 223, 
43 (1993); Tunneling in Complex Systems, edited by S. 
Tomsovic (World Scientific, Singapore, 1998) and refer- 
ences therein; A. Shudo and K. S. Ikeda, Phys. Rev. Lett. 
76, 4151 (1996); S. C. Creagh and N. D. Whelan, Phys. 
Rev. Lett. 77, 4975 (1996); E. Doron and S. D. Frischat, 



14 



Phys. Rev. Lett. 75, 3661 (1995); S. D. Frischat and E. 
Doron, Phys. Rev. E 57, 1421 (1998); J. Zakrzewski, 
D. Delande, and A. Buchleitner, Phys. Rev. E 57, 1458 
(1998). 

2" C. Eltschka and P. Schlagheck, Phys. Rev. Lett. 94, 014101 
(2005). 

3° J. U. Nockel and A. D. Stone, Nature 385, 45 (1997); W. 
K. Hcnsingor et. al, Nature 412, 52 (2001); D. A. Stock, 
W. H. Oskay, and M. G. Raizen, Science 293, 274 (2001); 
A. P. S. de Moura et. al, Phys. Rev. Lett. 88, 236804 
(2002). 

W. A. Lin and L. E. Ballentine, Phys. Rev. Lett. 65, 2927 
(1990); R. Utermann, T. Dittrich, and P. Hanggi, Phys. 
Rev. E 49, 273 (1994); D. Farrelly and J. A. Milligan, Phys. 
Rev. E 47, R2225 (1993); H. P. Breuor and M. Holthaus, 
J. Phys. Chem. 97, 12634 (1993); M. Gruebele, J. Phys. 
Condens. Matter 16, R1057 (2004). 

32 J. E. Baggot, Mol. Phys. 65, 739 (1988). 

33 S. Keshavamurthy and G. S. Ezra, J. Chem. Phys. 107, 156 
(1997); S. Keshavamurthy and G. S. Ezra, Chem. Phys. 
Lett. 259, 81 (1995). 

3-* M. Joycux, J. Chem. Phys. 109, 2111 (1998); M. Joyeux 
and D. Sugny, Can. J. Phys. 80, 1459 (2002). 

3^ M. Joyeux, S. C. Farantos, and R. Schinke, J. Phys. Chem. 
A 106, 5407 (2002); M. P. Jacobson and R. W. Field, 
J. Phys. Chem. A 104, 3073 (2000); H. Ishikawa, R. W. 
Field, S. C. Farantos, M. Joyeux, J. Koput, C. Bock, and 
R. Schinko, Annu. Rov. Phys. Chom. 50, 443 (1999). 

36 Z. M. Lu and M. E. KoUman, J. Chom. Phys. 107, 1 (1997). 

3^" J. Tennyson, N. F. Zobov, R. Williamson, O. L. Polyan- 
sky, and P. F. Bernath, J. Phys. Chem. Ref. Data 30, 735 
(2001). 

3® A. Bykov, O. Naumenko, L. Sinista, B. Voronin, J.-M. 
Flaud, C. Camy-Peyret, and R. Lanquetin, J. Mol. Spec. 
205, 1 (2001). 



33 S. Keshavamurthy, J. Chom. Phys. 119, 161 (2003). 

The cxpcriniontal value for the splitting of the (40)^4 pair 
is about 1.0 cm^ . J. Tonnysou, personal communication. 

A. J. Lichtenberg and M. A. Lieberman, Regular and 
Chaotic Dynamics, Springer- Verlag, New York, 1992. 
See, S. Keshavamurthy, J. Phys. Chem. A 105, 2668 (2001) 
for detailed references and one possible signature of non- 
linear resonances in parametric eigenvalue variations. 

''3 R. Pearman and M. Gruebele, J. Chem. Phys. 108, 6561 
(1998). 

B. V. Chirikov, Phys. Rep. 52, 263 (1979). 

A. Mouchet and D. Delande, Phys. Rev. E 67, 046216 
(2003). 

K. Husimi, Proc. Phys. Math. Soc. Japan 22, 264 (1940). 
See for example, A. Nitzan, Amm. Rev. Phys. Chem. 52, 
681 (2001) and references therein. 

M. Abu-Hilu and U. Peskin, J. Chem. Phys. 122, 021103 
(2005) and references therein. 

We note that it is perhaps better to think of Arnol'd dif- 
fusion in a slightly nonrestrictivc sense then that implied 
by Arnol'd in his seminal work, V. I. Arnol'd, Sov. Math.- 
Dokl. 5, 581 (1964). For critical reviews see: P. M. Cin- 
cotta, New Astr. Rev. 46, 13 (2002); P. Lochak in Hamilto- 
nian Systems with Three or More Degrees of Freedom, Ed. 

C. Simo, NATO ASI, Kluwer, Dordrecht, pp. 168, 1999; F. 
Vivaldi, Rev. Mod. Phys. 56, 737 (1984); G. Hallcr, Phys. 
Lett. A 200, 34 (1995). 

^° D. M. Leitner and P. G. Wolynes, Phys. Rev. Lett. 76, 216 
(1996); D. M. Leitner and P. G. Wolynes, Phys. Rev. Lett. 
79, 55 (1997); E. Tannenbaum, Ph.D. thesis, chapter 3, 
Harvard University, 2002. 

F. Grossmann, T. Dittrich, P. Jung, and P. Hanggi, Phys. 
Rev. Lett. 67, 516 (1991); J. Gong and P. Brumer, Annu. 
Rev. Phys. Chem. 56, advanced article, 2005. 



